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ABSTRACT 

We discuss the origin of the optical variations in the Narrow line Seyfert 1 galaxy 
NGC 4051 and present the results of a cross-correlation study using X-ray and optical 
light curves spanning more than 12 years. The emission is highly variable in all wave- 
bands, and the amplitude of the optical variations is found to be smaller than that of 
the X-rays, even after correcting for the contaminating host galaxy flux falling inside 
the photometric aperture. The optical power spectrum is best described by an unbro- 
ken power law model with slope a = 1.4^Q2 an d displays lower variability power than 
the 2-10 keV X-rays on all time-scales probed. We find the light curves to be signifi- 
cantly correlated at an optical delay of 1.2^ 'g days behind the X-rays. This time-scale 
is consistent with the light travel time to the optical emitting region of the accretion 
disc, suggesting that the optical variations are driven by X-ray reprocessing. We show, 
however, that a model whereby the optical variations arise from reprocessing by a flat 
accretion disc cannot account for all the optical variability. There is also a second sig- 
nificant peak in the cross-correlation function, at an optical delay of 39+g4 days. The 
lag is consistent with the dust sublimation radius in this source, suggesting that there 
is a measurable amount of optical flux coming from the dust torus. We discuss the 
origin of the additional optical flux in terms of reprocessing of X-rays and reflection 
of optical light by the dust. 
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1 INTRODUCTION 

The exact mechanism producing the X-ray emission from ac- 
tive galactic nuclei (AGN) is not well understood, although 
most models involve Comptonization of the UV/optical 
spectrum, through a tenuous, hot corona above the in- 
ner disc ((Haardt fe Maraschil Il99ll ). The fastest variabil- 
ity time-scales set a limit on the size of the emission re- 
gion, suggesting that the X-rays are emitted from the in- 
nermost regions of the system, but the exact geometry is 
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not known. The optical continuum is variable as well, and 
thought to be thermal emission from an optically thick ac- 
cretio n disc, surrounding the black hole l|Koratkar fc Blaea 
1999). The origin of the variability, as well as the connec- 
tion between the disc and the corona, is still unclear, how- 
ever. It has been suggested that the short time-scale opti- 
cal variations are the result of reprocessing of the highly 
variable X-rays by the disc dKrolik et al.l Il991f) . enhanc- 
ing the intrinsic thermal emission from the disc. Hence 
variations in the X-rays should be mirrored in the opti- 
cal light curves, smoothed by the large size of the disc 
and delayed by the light travel time between the source 
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and the disc. Although there are s ome sources that do 
show eviden c e of r eprocessing (e.g. IWanders et al.l 19971; 
Collier et aD Il998l: foknyanskii et all 120031 : ISergeev et all 
20051 ; ICackett et al.ll2007h . the results are not generally ap- 
plicabl e. Many sources show no simple rel ationship between 
bands (jNandra et aD l2000l : IciaskelfeoOSl) and at least on e 
source shows no correlation at all l|Maoz et al.l l2000. 2002). 
An increasing number of sources display optical variations 
of greater amplitude than t heir corresponding X-ray vari- 
atio ns, such as NGC 55 48 (lUttlev et all 120031) MR 2251- 
178 jArevalo et aill2008l) . Mrk 79 (|Breedt et al.ll2009l ) and 
NGC 3783 (|Arevalo etal]|2009t ). making it difficult to at- 
tribute all the optical variability to reprocessing. Energetics 
arguments also rule out reprocessing as the main driver of 
the variability in sources where th e variable opt ical luminos- 
ity exceeds the X-ray luminosity ([Gaskell 2 0071 ). 

Also important to consider, is variability intrinsic to 
the disc, such as instabilities propagating radially inward 
through the disc, affecting the mass accretion rate and lo- 
cal emission (e.g. lLvubarskiilll997n . In this case, optical 
variations are expected to precede similar variations in the 
X-rays. The time-scales involved are much longer, as vari- 
ations will propagate on the viscous time-scale, which is 
~ months-years at the optical emitting region in AGN. The 
long term (^months-years) optical variability has been ex- 
plained by composite models including a ccretion rate fluc- 
tuations propagating thr ough the disc (|Arevalo fe Uttlevl 



l200fj ; lArevalo et allboosl) and chan ges in the geometry of 



the system ( Breedt et alj 20091 ) . Reprocessed emission from 
the broad line region may also contribute to the o ptical 
variatio ns ( Arevalo et al.l [20091 ; iKorista fc G oad 2001), and 
iGaskelll l|2007l ) suggested that there might be a contribution 
to the optical emission from reprocessing material much fur- 
ther out, such as the molecular torus. 

By undertaking simultaneous multi-wavelength moni- 
toring and studying the correlation between the emission at 
different wavelengths, we can determine the dominant pro- 
cess underlying the optical variations and gain some physical 
insight into the connection between the emission regions. 

lUttlev et al.l (|2003l ) proposed that the dominant pro- 
cess driving the observed optical variability, was depen- 
dent on the mass and accretion rate of the black hole, as 
these parameters define the temperature profile of the ac- 
cretion disc. Within the st andard accretion disc model (e.g. 
IShakura fe Sunvaevl Il973l ). AGN with smaller black hole 
mass and/or higher accretion rate, will have have hotter 
discs than systems of high mass/low accretion rate. Scaled 
in terms of the gravitational radius, the X-ray and optical 
emitting regions are further apart in the low mass systems, 
so the viscous time-scale in the optical emitting region is 
very long compared to that in the inner disc. We may there- 
fore expect reprocessing to be the main contributor to the 
fast optical variations, and accretion rate variations to occur 
on time-scales of several tens of years. 

NGC 4051 is a low-luminosity narrow-line Seyfert 1 
galaxy, known for its extreme X-ray variability (e.g. 
iMcHardv et al.ll2003 ). The relationship between the X-ray 
and UV/optical variability in this galaxy has been the sub- 
ject of many papers, reporting varied and conflict ing re- 
sults. Over the course of 3 nights iDone et al.l l|l990l ) found 
the X-rays to vary by a factor of 2, but no corresponding 
variations in the optical or infrared were seen. iKlimek et al.l 



l|2004l ) also studied the short time-scale optical variability of 
this galaxy. They detect ed microvariability o n one of the five 
nights it was observed. lUttlev et all (|2000h found a strong 
correlation between the X-ray and EUV varia tions, and mea- 
sured them to be simultaneous to within 1 ks. lPeterson et all 
(2000) concluded that the X-ray and optical emission was 
correlated on time scales of months to weeks, but by ap- 
plying a 30 day smoothing boxcar to the light curves be- 
fore calculating the cross-correlation, a measurement of any 
lag shorter than this could not be obtained. During a pe- 
riod of very low X-ray flux they found that the broad He 
11 A4686 line greatly decreased in intensity as well, show- 
ing no noticeable variation during this period, but that the 
optical continuum and H/3 line were only slightly fainter 
and continued to vary. They interpreted this is as the inner 
disc turning into an advection-dominated flow, greatly re- 
ducing the short wavelength (UV/blue) emission from this 
part of the disc. A 3-mont h simultaneous X-ray and opti- 
cal monitoring program bv lShemmer et all (|2003l ) revealed 
significant evidence of a correlation close to zero lag. They 
found the optical emission in this case to lead the X-ray 
emission by 2.4 days, and note that the variability of the 
optical light curve was much less than tha t of the X-rays. 
On shorter time-scales, iMason et all (|2002h find the UV to 
lag the X-rays, which they interpret as reprocessing by a 
ring of material at a distan ce of 0.14 light-days. Thi s result 
has also been confirmed bv lSmith fc Vaughanl l|2007h . using 
the s ame data. Reverbera t ion mapping studies using the H/3 
line ilPeterson et al.ll200ol ; iKaspi et aDl2000l ; iPeterson et~atl 
|2004| ) place the broad line region in this AGN at a distance 
of ~6 light-days from the central ionising source. 

ft is clear that a combination of processes can con- 
tribute to the observed optical variability and that differ- 
ent processes can drive the variability on different time- 
scales. In this paper, we combine the observations from 
4 different telescopes with some archival data, to revisit 
the topic of X-ray and optical variability in NGC 4051. 
The light curves together span more than 12 years, al- 
lowing us to study the variability on much longer time- 
scales than has been done before in this source. This pa- 
per is the fourth in a series describing our findings from 
a sample of Seyfert galaxies monitored at X-ray and opti- 
cal wavelengths. Previously we have reported on our results 
on M R2251-178 (jArevalo et aDl2008l), Mrk 79 (|Breedt et all 
120091) and NGC 3783 |Arevalo et alj|2009l ) and here we dis- 
cuss our findings on N GC 4051, the lowest mass (M = 
(1.91 ± 0.78) x 10 6 M^ . [Peterson et al.ll2004l ). lowest lumi- 
nosity (Lboi < 10 44 erg s _1 ) galaxy in our sample. Through- 
out we ad opt the avera ge Tully-Fisher distance estimate of 
15.2 Mpc (|RusseUll2003l ) to this source. 

In Section [2] we present the multicolour optical and 
X-ray data. We show the optical power spectrum in Section[3] 
and in SectionUwe present the results of a cross-correlation 
analysis of the light curves. We discuss the results in terms 
of a simple reprocessing model in Section [5] and show that 
there is an additional component to the optical variability. 
We show that these additional variations are delayed with 
respect to the X-rays by a time-scale consistent with the 
light travel time to the torus, so in Section [6] we investigate 
the possibility of optical emission originating from the torus, 
in detail. Our findings are summarised in Section [7] 
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2 DATA 

2.1 V band optical observations 

The V band optical light curve presented in this paper 
is a combined light curve from four telescopes: the Liver- 
pool Telescope (LT) at the Observatorio del Roque de Los 
Muchachos of t he Instituto de As trofisica de Canarias, La 
Palma, SpairQ l|Steele et al.l [2004). the Faulkes Telescope 
North (FT) at Haleakala, Maui, Hawaii, the MAGNUM 
(Multicolor Active Galactic Nuclei Monitoring) telescope 
jYoshii et al.ll2003l ) at the University of Hawaii's Haleakala 
Observatory, Maui, Hawaii, and the 0.7-m telescope of the 
Crimean Astrophysical Observatory (CrAO) in t he Crimea, 
Ukraine i|Sergeev et al.ll2005l ; ICnesnok et al.ll2009l ). The data 
reduction procedur es and light curve c onstruction are de- 
scribed in detail in IBreedt et ail (|2009l) , so we only give a 
brief summary below. 

We observed NGC 4051 through a Bessell-U band filter 
from 2004 December 9 to 2008 December 11 with the Liver- 
pool Telescope and from 2007 February 5 to 2007 August 10 
with the Faulkes Telescope North. The two telescopes are 
identical in design. Observations in other filter bands, on 
both telescopes, started in early 2007 and continued until 
2007 August on the FT and 2008 July on the LT. The data 
reduction and photometry was done using standard tech- 
niques in iraf. We performed aperture photometry on the 
individual images, using an aperture of 12" diameter, cen- 
tred on the nucleus of the galaxy. The relative fluxes were 
calibrated with respect to four comparison stars in the field, 
and the fi nal flux calculated u sing t he comparison star mag- 
nitudes of IDoroshenko et al.l (120051 ). 

The MAGNUM telescope observed NGC 4051 in the V 
band from 2001 January 8 to 2007 July 30. Data from the 
start of this progr am until MJD 52819 (20 03 June 29) have 
been published by ISuganuma et al.l (|2006l ). The rest of the 
light curve was previously unpublished. 

Observations at the CrAO were carried out between 
2001 December 27 and 2006 August 18. The filter set on the 
Crimean telescope is non-standard, but their B,V,R and / 
filters are close to the Johnson/Bessell fi lters and their Rl 
filter is similar to the Cousins I filter (see IDoroshenko et al.l 
120051 ). 

The light curves from the different telescopes were made 
using different aperture sizes, so each is affected by a differ- 
ent amount of host galaxy contamination, resulting in dif- 
ferent average flux levels. To account for this, as well as for 
the slight differences in the filters and any other small sys- 
tematic differences in the calibration, we use the overlapping 
parts of the light curves to match them in a least squares 
fashion. We extracted pairs of measurements, taken within 
one day of each other, from overlapping light curves, and 
performed a least-squares fit of the equation Fi = aF-z + b 
to the flux pairs. F\ and F2 axe the flux measurements from 
the two datasets and a and b the parameters of the fit. The 
points were weighted for the fit, by the inverse of the time 
difference between Fi and F2, so that measurements made 
closer together in time have a greater influence on the fit. 



1 The Liverpool Telescope is operated by Liverpool John Moores 
University with financial support from the UK Science and Tech- 
nology Facilities Council. 



The best-fit parameters from each pair of overlapping light 
curves were then used to scale that dataset to the 15" diam- 
eter aperture V band light curve from the CrAO telescope. 

The V band light curve used in the analysis, is the final 
combined light curve, corrected for Galactic reddening and 
host galaxy flux inside the aperture (see Section 12. 5p . It is 
shown in Figure [1] normalised to its mean flux, for direct 
comparison with the X-ray variations. Full sampling char- 
acteristics of this light curve may be found in Table [T] 

2.2 AGN Watch continuum light curve 

NGC 4051 was also observed spectroscopically between 1996 
January 12 and 1998 July 2 8 as part of th e International 
AGN Watch program (e.g. IPetersonI Il999l ). Light curves 
from this program are publicly availabl(Q and we include 
in our analysis t he 5100A continuum light curve origi- 
nally presented bv lPeterson et all (|2000l ). We subtracted the 
host galaxy contribution inside the 5.0" x 7.5 " spectroscopic 
aperture, as measured bv lBentz et ail l|2009l ). and converted 
the light curve to mjy units, to compare with our V band 
light curve. Note that in this light curve the continuum flux 
is measured from a line-free region around 5100A of the op- 
tical spectra, while the broad band observations will contain 
the line emission that fall inside the V band bandpass, such 
as [Om]AA4959,5007 and H/3A4861. We discuss the contri- 
bution of the lines to the flux in Section ^. II The light curve 
is plotted in Figure [TJ together with the X-ray and V band 
light curves. 

2.3 Multicolour optical observations 

In addition to the V band observations described above, the 
CrAO also observed NGC 4051 in the B, R, Rl and I filter 
bands during the same period. The light curves are con- 
structed from aperture photometry through a 15" diameter 
aperture centred on the galaxy nucleus, relative to a com- 
paris on star in the field. We use the comparison star fluxes 
from lCackett et ail (2007) to convert the light curves to flux 
units in these bands. The data from the start of this pro- 
gram, unt il MJD 53022 (2004 January 18), have been pub- 
lished by ISergeev et ail (|2005l ) and were also analysed by 
ICackett et al.l (|2007l ). Data after this date were previously 
unpublished. 

We supplement these light curves with Sloan Digital 
Sky Survey (SDSS)-u band observations made by the Liver- 
pool Telescope between 2004 December 9 and 2008 July 27. 
The data from short monitoring programs in B and R on the 
Liverpool and Faulkes telescopes (2007 February 5 - 2008 
July 27) were combined to extend the CrAO light curves. 
Unfortunately there is no overlap between the LT/FT and 
the CrAO light curves in these bands, so it was not possi- 
ble to combine them as described in section 12.11 To ensure 
the best possible match between the light curves, the same 
comparison stars were used in the flux calibration of the two 
datasets, and the appropriate aperture corrections were ap- 
plied to the LT/FT light curves. However, it is possible that 
small differences in the calibration, or between the filters on 
the different telescopes, may cause a small offset between 

2 http:/ /www. astronomy.ohio-state.edu/~agnwatch/ 
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Figure 1. 3 — 10 keV X-ray (top panel in each part) and optical (bottom panel in each part) light curves of NGC 4051, corrected for 
host galaxy flux inside the aperture. Each light curve is normalised to its mean flux for direct comparison of the variability amplitudes. 
Note that the vertical scales for the optical and X-ray light curves are different — the amplitude of the optical variations is smaller than 
that of the X-rays. The X-ray light curve is smoothed with a 4 point (~1 week) running average (for plotting purposes only), to highlight 
the variations on time-scales similar to that of the optical light curve. The first part of the optical light curve, plotted in blue triangles, 
is the 5100A continuum light curve from the AGN Watch nroeram. The rest of the light curve, nlotted in green dots, is the combined 
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Figure 2. X-ray, u, B, V, R, Rl and I band light curves of NGC 4051. The X-ray light curve is in units of 10 — 11 erg s _1 cm -2 and 
the optical light curves are in mjy. The fluxes have been corrected for the small amount of Galactic extinction towards the source, but 
host galaxy flux has not been subtracted. The dotted lines in the B and R band light curves indicate the positions where there might be 
a small offset between the light curves on either side of the line, as there was no overlap between the data from the different telescopes. 
This does not affect our cross-correlation results however (see text). 



Table 1. Observational characteristics of all light curves 



Light curve 


Telescope 


Observation dates 


Total length 


Number 


Ai mcan 


At mcd ia„ 


Mean flux* 






(MJD) 


(d) 


of epochs 


(d) 


(d) 


(mjy) 


X-ray 


RXTE 


50196 - 54811 


4615 


1811 


2.55 


1.96 


1.93 ±0.95* 


u 


LT 


53348 - 54675 


1327 


128 


10.4 


4.01 


5.41 ±0.93 


B 


combined 


52271 - 54675 


2404 


386 


6.24 


1.96 


11.5 ±0.61 


5 100 A 


AGN Watch 


50095 - 51023 


928 


126 


7.42 


3.15 


11.1 ±0.93 


V 


combined 


51917 - 54811 


2894 


629 


4.61 


1.94 


20.1 ±0.64 


R 


combined 


52271 - 54675 


2404 


384 


6.28 


1.97 


32.3 ±0.91 


Kit 


CrAO 


52271 - 53965 


1694 


265 


6.42 


1.98 


38.9 ± 1.02 


/ 


CrAO 


52271 - 53965 


1694 


261 


6.52 


1.98 


41.3 ± 1.02 



* Including host galaxy flux 

^ In units of 10 — 11 erg s _1 cm -2 

t The Kl filter corresponds approximately to the Cousins I filter 
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the CrAO and LT/FT light curves. The positions in the 
light curves that are affected by this, are indicated by dotted 
lines in Figure [2] Note however, that due to the "segmen- 
tation" technique we use to calculate the cross-correlation 
(Section a possible small offset like this does not affect 
our results. The light curves are shown, together with the 
corresponding segment of the X-ray light curve, in Figure [2] 
The flux uncertainty on each point, based upon the dif- 
ferential photometry error and the relative calibration of the 
light curves, is approximately 3%. The uncertainty in the 
comparison star fluxes introduce a 1 — 4% error in average 
flux of the different bands, but note that this error does not 
increase the scatter in the light curve, so it does not affect 
the cross-correlation results below. Sampling statistics for 
each light curve are shown in Table [T] 

2.4 X-ray observations 

NGC 4051 is part of our long term AGN X-ray monitor- 
ing program using RXTE, which started on 1996 April 23. 
Here we include the observations up to 2008 December 11, to 
overlap with the optical light curve. The observational fre- 
quency and variability characteristics of this light curve are 
summarised in Tables Q] and [5] The observations are in the 
form of approximately 1 ks snapshots using the Proportional 
Counter Array (PCA) on board RXTE. The PCA consists 
of five identical proportional counter units (PCUs) but as 
some of them were regularly switched off during our ob- 
servations, we only extracted data from PCU2 to construct 
the X-ray light curve. For maximum signal-to-noise, we used 
only the data from the top layer. Using FTOOLS v. 6. 4, we ap- 
plied standard extraction techniques and re-reduced all the 
available data, so that the latest RXTE background models 
may be applied throughout the light curve. We rejected data 
obtained within 10 minutes of passing through the South At- 
lantic Anomaly (SAA), data gathered less than 10° above 
the limb of the Earth or with pointing offset greater than 
0.02° from the source, and data for which the electron cur- 
rent was greater than 10%. To obtain the flux measurement, 
we fit a power law model to the spectra using xspec and in- 
tegrate the flux in the range 3-10 keV. The error on the 
flux measurements were calculated from the counts in each 
observed spectrum, weighted by the response function. 

We show the normalised 3-10 keV X-ray light curve 
in Figure [T] smoothed with a 4-point (~1 week) running 
average. This is done to highlight fluctuations on similar 
time-scales as those in the optical light curve, but note 
that, throughout this paper, the unsmoothed, unbinned light 
curve was used in the analysis. 

2.5 Galaxy flux in the aperture 

The photometry aperture contains the nuclear emission as 
well as a constant contribution from stars in the host galaxy. 
Although this constant offset does not affect the results of 
the cross-correlation analysis, it reduces the amplitude of 
the observed variations. Hence, to investigate the variability 
of the nucleus, this component must be subtracted from the 
light curve. 

In order to estimate the amount of host galaxy flux con- 
tained in th e 15" aperture, w e decomposed the images us- 
ing GALFIT |Peng et al.l l2002). by fitting a psf (point spread 



function - using stars in the field of the galaxy), a de Vau- 
couleur profile and an exponential disc to each image. The 
decomposition of individual images was poor, so we stacked 
together 12 images taken under good seeing conditions, to 
improve the signal-to-noise. This resulted in much better, 
stable fits of the different components. The variable nucleus 
was taken account of by selecting only images from which 
we measured approximately the same flux through the aper- 
ture. We subtracted the best fit galaxy profile in each band 
from the images and measured the flux of the residual psf- 
like image. The difference between this psf- like image and 
the flux measured through a 15" aperture on the original 
image, was taken as the galaxy contribution f g . We list these 
fluxes, corrected for foreground Galactic extinction, in Ta- 
bled 

We did not have Rl and / band images available to mea- 
sure the galaxy contribution by decomposition as above, so 
the galaxy contribution in these bands were estimated by 
extrapolating from the galaxy contribution measured in the 
SDSS-i band. The fluxes in Table [2] arc in good agreement 
with those found bv lSalvati et al.l (1993) and also agree well 
wi th the galaxy-subtracte d V band light curve published 
by ISuganuma et al.l (|2006l ). We show in Figure [3] the mea- 
sured galaxy fluxes, converted to F\ units and corrected 
for the small amount of foreground Galactic extinction to- 
wards the source. It is shown together with a template spec- 
trum of the centr al regions of a normal spiral g alaxy simi- 
lar to NGC 4051 (ISantos et al.ll2002l ; lBicalll988T l. The tem- 
plate spectrum was constructed from a group of 30 galaxies, 
mostly of morphological type Sb, but also including Sa and 
Sc galaxies. NGC 4051 is classified as an Sbc galaxy. The 
template is corrected for Galactic extinction and normalised 
at 5870A. Our galaxy fluxes agree well with the template 
spectrum, indicating that our f g measurements are good es- 
timates of the host galaxy flux inside the aperture. Flux-flux 
plots between galaxy-subtracted light curves in different en- 
ergy bands suggest that the error on f g is between 10 and 
15%. 



2.6 Fractional variability 

The fractional variability F val describes the flux variance of 
the light curve as a fraction of the mean flux (/), and pro- 
vides a way of comparing the amount of intrinsic variability 
in each wavelength band. Mathematically, 



Here a is the flux variance 

JV 

^£(/*-</» 2 



(i) 



(2) 



and (e ) is the mean of the squared measurement errors 

N 



The error on F var is given by 



£ (F var ) - 



\ 




(3) 



V 777 (4) 
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Wavelength [A] 

Figure 3. Measured galaxy flux inside the aperture (plotted in 
dots), compared to a template spectrum of the central regions 
of a normal spiral galaxy similar to NGC 4051. The last two 
points, plotted in squares, are the flux estimates for the Rl and 
/ bands. These points were extrapolated from the SDSS-i flux 
measurement (last dot), as we did not have images available in 
these bands. 

ijVaughan et alj l 2003) . This error describes only the error on 
F val . due to observational noise. In our case, there is also an 
uncertainty associated with the galaxy subtraction. This un- 
certainty does not affect the scatter of the points in the light 
curve (since we are subtracting the same galaxy contribu- 
tion to each flux measurement), so should not be included in 
the average error in Eqn.[3] It should be accounted for in the 
error on F v!a however, so we also include a third error term 
in the square root of Eqn. H]of the form (F v „ A(/)/{/) ) 2 , 
where A(/) is the error in the galaxy flux contribution. 

The fractional variability, as well as the ratio of the 
maximum to minimum flux in each light curve, is given in 
Table [2] We also list the variability characteristics of the 
5100A continuum light curve togeth e r with the galaxy con- 
tribution measured by iBentz et al.l (|2009l ). but point out 
that this light curve was made using a smaller aperture than 
our light curves and that the galaxy contribution was mea- 
sured from much higher resolution images than we used here. 
There may therefore be a difference in the amount of host 
galaxy flux subtracted, so the fractional variability in this 
band should not be compared directly with the other light 
curves. The X-ray and optical light curve s are all highly 
variab le. In agreement with the finding of She mmer et al.l 
l|2003h that the optical continuum variations are less than 
the X-rays on time-scales of months, we find that the same 
is true on time-scales of years. The optical continuum (all 
light curves, including the 5100A light curve) is less vari- 
able than the X-rays, even after correcting for host galaxy 
starlight in the aperture. 

In the X-ray reprocessing scenario, the fast variations 
of the X-ray light curve are smoothed out by the disc, so a 
smaller fractional variation of the optical light curves is ex- 
pected. Our results are broadly consistent with this expecta- 
tion, showing a decrease in F var towards longer wavelengths. 
However, due to the uncertainty in the galaxy subtraction, 
the error on F mi is large, and the decreasing trend is not 
strong. The observation that the variability of the optical 



Table 2. Estimated host galaxy flux inside the 15" aperture and 
variability characteristics of all light curves, after subtracting the 
host galaxy contribution. 



Light curve 


fg (mJy) 


F var (%) 




X-ray 




48.8 ±0.8 


89.0 ±26.5 


u 


1.8 


25.6 ±1.8 


3.92 ± 1.36 


B 


6.3 


11.7 ±1.5 


1.92 ± 0.56 


V 


16.0 


15.4 ±3.1 


2.48 ± 1.64 


R 


27.2 


16.9 ±3.8 


2.46 ± 1.65 


Rl 


31.5 


13.8 ±3.3 


2.02 ± 0.92 


I 


32.7 


11.8 ±3.0 


1.84 ±0.72 


5100A 


7.2* 


23.9 ±1.5 


2.73 ± 0.15 



jBentz et alj l l2009h . in the 5.0" x7. 5" spectroscopic aperture. 



bands is less than that of the X-rays is nevertheless robust. 
We note that the variability of the B band light curve is 
lower than the other optical bands and does not follow the 
decreasing trend. The mean flux of the galaxy-subtracted B 
band light curve is also slightly higher than that of the V 
band, which is contrary to expectation, assuming the emis- 
sion in both bands originate in an accretion disc. This sug- 
gests that the low fractional variability in the B band may be 
the result of underestimating the galaxy contribution in this 
band. Increasing the galaxy contribution by 20%, increases 
F var to 15.5 ± 1.3%. The larger error may be attributed to 
the fact that the galaxy spectrum rises steeply in the B 
band (see Fig. [3} , and is likely to be very different from the 
spectrum of the comparison star used in the differential pho- 
tometry. The lower variability of the optical with respect to 
the X-ray variations, suggests that reprocessing of X-rays 
may account for the optical variability in this source. We 
return to this possibility in Section [5] 



3 THE OPTICAL POWER SPECTRUM 

Power spectra are a commonly used technique to investigate 
the aperiodic variability displayed by AGN and X-ray bina- 
ries. It describes the variability power present in the light 
curve (mean squared amplitude) as a function of frequency. 
Generally, the X-ray Power Spectral Density (PSD) of these 
systems can be described by a power law P(y) cx v~ a of 
slope a ~ 1 up to a bend or break frequ ency vb, where 
the power law slope s teepens to a ~ 2 (e.g. iMcHardv et al.l 
l2004l : ISummonsll2008h . The corresponding break time-scale 
scales approximately linearly with the black hole mass, and 
the scatter in this relationship is largely accounted for by 
consideri ng the difference in accretion rate between these 
systems (|McHardv et al.ll2006l ). The break frequency is thus 
characteristic of the system and is thought to represent a 

time-scale associated with the inner edge of th e disc. 

W e us e the Mo n te Ca rlo technique of lUttlev et al.l 
(2002) and S ummonsl (|200ct ) to estimate the parameters of 
the power spectrum underlying the optical variations. The 
method subtracts the mean of each of the component light 
curves and then calculates the discrete Fourier Transform of 
each part. The resulting PSD is binned logarithmically in 
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bins of width 1.5v, with v the frequency at the start of the 
bin. 

We do not include the AGN Watch data here, as the 
slightly different wavelength will affect the normalisation of 
the power and hence the slope we measure. To extend the 
frequency r ange, we also includ e the V band light curve 
obtained bv lKlimek et al.l (|2004T ) as part of their microvari- 
ability study. We calculate the variability power on each 
of the five nights it was observed and bin the resulting 
PSDs together to give the power in the frequency range 
6 x 1(T 5 - 5 x 10~ 4 Hz. The low frequency part of the PSD 
is calculated from the whole of the V band light curve and 
the middle range from the most intensively sampled part of 
the light curve from MJD 54181 to 54252. 

We start by fitting an unbroken power law 



P(v) = Au~ 



(5) 



to the power spectrum, allowing the slope a and the normal- 
isation A to vary. The best fit, shown in Figure f?] is found 
to have a = L4fto"a, with an acceptance probability of 79.5 
per cent. We also include the shape of the X-ray power spec- 
trum, plotted as a dashed line, for comparison. The X-ray 
power spectrum has a slope at = 1-lio' \ at low frequen- 
cies, bending to a slope of an = 2.5^°'° and a frequency 
v B = 5-llte x 10 ~ 4 Hz (|Summonj|200a Summons et al., 
in preparation). 

The unbroken optical PSD slope is intermediate in value 
between the low and high frequency X-ray slopes. If the op- 
tical PSD has a similar shape to the X-ray PSD, this sug- 
gests that a break is contained within the frequency range 
spanned by the optical data. Motivated by the shape of the 
X-ray PSD, we also fit a bending power law, of the form 



P{y) 



Av~ aL 



1 + (v /v B ) aH ~ aL 



(6) 



to the optical PSD. For direct comparison with the X-ray 
PSD, we fixed the low frequency slope to ai = 1.1 and 
allowed the high frequency slope an and bend frequency 
vb to vary. We find the best fitting parameters to be an = 



1.5? 



and bend frequency vb = 1.4_, x 10 Hz. The 



errors are the values for which the acceptance probability 
drops below 10 per cent and the asterisk is used to indicate 
an unbounded error. The acceptance probability for this fit 
is slightly better than for the unbroken case, at 90.2 per cent, 
but the break frequency is unconstrained by the frequency 
space sampled by the data. Allowing all the parameters to 
simultaneously vary over a large range (i.e. not fixing the 
value of the low frequency slope), results in the same best- 
fit values, within the quoted 90% errorbars. 

We show in Figure [5] a colour density plot of the ac- 
ceptance probability as a function of the parameter space 
searched, while fixing ql = 1.1. It shows that for a.u = 2.5, 
as measured from the X-ray spectrum, the best-fit bend fre- 



quency is at vb 



■■ 4.4t* 2 x 10" 



Hz. The fit probability for 



these parameters is only 74.7 per cent, however, so it can be 
confidently rejected in favour of the unbroken model. It is 
clear from this figure that the best-fit bend frequency in the 
single-bend model is towards the lowest frequencies (longest 
time-scales) probed. A single-bend PSD cannot be ruled out 
for this system, but it would require further monitoring to 
confirm the existence of a bend at low frequencies. We there- 
fore conclude that, with the current data, the optical power 




Figure 4. Unbroken power law model (solid line) fitted to the 
unfolded optical PSD (points with error bars). The slope of the 
best-fit model is a = 1.4. The shape of the X-ray power spectrum 
jSummonsll20u8h is shown in a dashed line for comparison. 



log v B (Hz) 




Figure 5. Acceptance probability as a function of oljj and vb 
for the parameter space searched, was fixed at 1.1, for direct 
comparison with X-ray PSD. The best fit parameters yield an 
acceptance probability of 0.902. The lines superimposed on the 
colour plot are the 67% and 90% confidence contours. 



spectrum is best described by an unbroken steep power law 
model with slope a = 1.4. 



4 CROSS-CORRELATION ANALYSIS 
4.1 X-ray — optical correlation 

In order to test whether there is a measurable delay of one 
light curve with respect to the other, we calculate their cross- 
co rrelation function ( CCF) , usin g the interpolation metho d 
of iGaskell fc Sparkd (I1986T ) and I White fc Peterson! (|l994h . 
We implement the method in two ways: in the first case, we 
linearly interpolate between adjacent points and resample 
both light curves to obtain equally spaced points for calcu- 
lating the correlation. We will refer to this method as the 
"equally sampled" (EQ) method. We note however that for 
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most of the monitoring time, the X-ray light curve is bet- 
ter sampled than the optical light curve. To take advantage 
of this and maximise the use of real, rather than interpo- 
lated data, we use a second implementation of the method, 
where we interpolate only the optical light curve and use 
the observed fluxes and observation times of the X-ray light 
curve in the calculation. We will designate this the "single 
interpolation" (SI) method. The two methods give consis- 
tent results throughout, but for completeness we list the 
cross-correlation results from both methods in Table [3] 

To take account of the gaps in the optical light curve 
when the galaxy is not observable from ground-based ob- 
servatories, we divide the light curve into segments, defined 
by these gaps. We calculate the correlation between each 
segment and its corresponding segment of the X-ray light 
curve and average the resulting CCFs, weighting them by 
the length (in time) of each segment. The final CCF, calcu- 
lated in this way, is shown in Figure [S] It shows a clear peak 
at 1.21q'3 days, and another at 38.9lg 4 days. The errors 
are calculated using a simple bootstrap method. We select 
at random two-thirds of the points in the X-ray light curve, 
each along with its closest optical data point, recalculate the 
CCF and then measure the peak and centroid. We performed 
1000 such selections to yield a distribution of peak and cen- 
troid values. The error range reported contains 68% of the 
measured values about the median. Here, and throughout 
this paper, a positive lag indicates the longer wavelength 
band lagging behind the shorter wavelength band. 

The segmentation technique also allows us to take ac- 
count of the small difference in the flux level of the V band 
and 5100A light curves due to different aperture sizes and 
a possible difference in the galaxy subtraction. The mean of 
each segment is calculated individually and subtracted from 
the segment before the correlation is calculated. The seg- 
mentation has the effect of removing long time-scale power 
in the light curve and allows a more accurate measure of 
the short time-scale lag between the light curves. Any small 
differences in the relative calibration of the light curves are 
therefore eliminated from the calculation. 

We perform Monte Carlo simulations to test the sig- 
nificance of these lags in the following way: we generate 
1000 random red noise light curves with the same statis- 
tical properties a s the observed X-ray lig ht curve, based 
on the method of iTimmer fc Koenid (Il995l) and the X-ray 
power spectrum parameters of NGC 4051 (|Summonsll2008l . 
see also this paper, section [3}. Each of these random, mi- 
correlated light curves is then sampled in the same way as 
the observed X-ray light curve and cross-correlated with the 
real, observed optical light curve, using the same segmen- 
tation technique as described above. The resulting CCF is 
recorded in each case and compared with the CCF of the ob- 
served light curves. In this way we can assess the probability 
of finding spurious correlations, i.e. chance correlations due 
to the red noise character of the light curves, that are higher 
than that of the real light curves. We overlay the mean, 95% 
and 99% levels of these simulations on top of the calculated 
CCF, in Figure [5] Both peaks reach higher than the 99% 
line, showing that the correlations are significant at greater 
than 99% confidence. 

To illustrate the origin of the two lags in the light curves, 
we show in Figure[7]the optical light curve shifted by 1 (top) 
and by 39 days (bottom), in comparison to the X-ray light 
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Figure 6. CCF between the X-ray and long term optical light 
curves, shown in Figure[TJ The horizontal dotted, dot-dashed and 
dashed lines are the mean, 95% and 99% confidence levels, cal- 
culated from 1000 Monte Carlo simulations. The vertical dotted 
line indicates the zero lag position. Both the peaks are significant 
at greater than 99% significance. 
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Figure 7. The optical V band light curve shifted back by 1 d 
(top) and 39 d (bottom), compared to the X-ray light curve (mid- 
dle). The X-ray light curve was binned in 5 day bins for clarity. 
The short time-scale correspondence between X-ray and the opti- 
cal light curve at the top is clear. The vertical lines indicate some 
of the variations which align better when the optical light curve 
is shifted back by 39 days, rather than by the main 1 day lag. 



curve (middle). The X-ray light curve was binned in 5-day 
bins for clarity. Most of the X-ray variations on time-scales 
of days to weeks show corresponding variations in the optical 
light curve shifted by 1 day. However, a few of the minima 
and maxima align better with the 39-day shifted version of 
the optical light curve. We point out that since both lags are 
short compared to the total duration of the light curve, only 
the largest amplitude variations can be clearly seen in Fig- 
ure [7] The most obvious examples of the delayed variations 
are marked by dotted lines in that figure. 
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4.2 Lags between optical bands 

It is clear from the light curves in Figure [2]that a strong cor- 
relation exists between the optical bands. Opt ical interband 
lags i n this source have b een m easured by ISergeev et al.l 
(2005) and ICackett et all |2007l ). confirming the apparent 
correlation. We consider here the relationship between the 
X-rays and each of the optical bands. 

Using the same segmentation technique as before, we 
calculate the correlation between the X-ray light curve and 
each of the different optical bands. The CCFs are shown 
in Figure [S] and we summarise in Table [3] the value of the 
peak and centroid of the CCF, calculated with both methods 
of interpolation. The centroid is calculated as the weighted 
average of the CCF points around the peak, at a value 
^ 85% of the peak correlation coefficient. As before, a pos- 
itive value of the lag indicates that the variations in the 
optical band are lagging behind the X-ray variations. All 
light curves, except the u band, for which the sampling is 
much sparser than the other light curves, show a correlation 
significant at greater than 99% confidence. The second peak 
near ~ 40 days is also seen, though not always at a signifi- 
cant level. We note however, that the strength of the second 
peak, relative to the peak at ~ days, remains approxi- 
mately constant. The V band light curve used in this part 
of the analysis, consists of only those observations which 
overlap in time with the other colour light curves. From 
this segment of the data we measure a slightly shorter lag 
than from the full light curve, but the results are consistent 
within the la errors. We also note that the lags measured 
from the central CCF peak are well within the ~ 6 light 
day estimate of the location of the H/3 broad line region 
(BLR) in this source (|Peterson et alfcOOfj ; iKaspi et al1l2000l ; 
iPeterson et~ail I2004T I . indicating that most of the variable 
optical emission originates in a structure smaller than the 
BLR. 

In order to assess the influence of the delayed line emis- 
sion on the measured lag, we first estimate the line contri- 
bution to the flux measured through the broad band filters. 
Using an average optical spectrum from the AGN Watch 
program, convolved with the filter response, we find a 3% 
line contribution to the V band flux and a 4% contribu- 
tion to the B band flux. We then generate correlated X-ray 
and optical light curves, using the method described in Sec- 
tion UT] to represent the continuum emission. The line flux 
is represented as a smoothed version of the X-ray light curve, 
delayed by 6 days, as would be appropriate for a thin spheri- 
cal shell of radius 6 light-days. The delayed line flux is added 
to the simulated, undelayed optical light curve in the ratios 
measured from the spectrum. 1000 such simulations yield 
a peak lag of r pcak = 0.0 ±0.1 days and a centroid lag of 
T cont = 0.1 ± 0.6 days, which is well within the errors on the 
lag meaured in the real light curve. We therefore conclude 
that the line emission contributes too small a fraction of the 
total V band emission to noticably affect the measured lag. 



5 REPROCESSING 

The short, but measurable delays between the X-rays and 
the optical bands, as well as the small amplitude of the op- 
tical variations, suggest that reprocessing of X-rays plays 
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Figure 8. CCF between the X-ray and optical bands as labeled. 
The vertical dotted line indicates the zero lag position. For sim- 
plicity, we only show the 95% (dot-dashed line) and 99% (dashed 
line) confidence levels. 



Table 3. Cross-correlation lags (in days). The centroid is the 
weighted average of the CCF values ^ 85% of the peak correlation 
coefficient. 



CCF 
X-ray vs. 


Peak 
(EQ) 


Centroid 
(EQ) 


Peak 
(SI) 


Centroid 
(SI) 


Optical 11 












first peak 


2.0 


fl.2 
-0.6 


2 4+ - 9 

z '^-1.6 


1 2 +1 -° 


1 8+ 1 ' 3 
i - S -2.0 


second peak 


35.2 


+6.6 
-5.1 


39.4+J; 2 


38 - 9 ±8.I 


38.7±^;g 


u 


Q.3 


fl.5 
-1.3 


5+ 1 ' 4 
u - a -i.6 




l +2 '° 


B 


2.3 


fl.5 
-0.7 


2 6+ - 5 


i =+0.8 
1,o -0.6 


2.0±°-* 


V 


1 1 +o.s 

J - L -0.3 


9 9 +0.3 
Z - Z -0A 


6+ - 4 


i q+0.8 
1 " 3 -1.0 


Jrl 


3.1 


fl.O 
-1.6 


2 8+ 06 


1 

i '°-0.6 


2.o±;-; 


Rl 


4 8+2-1 


n q + 1.3 

Z - 3 -1.7 


1 9+ 31 




1 


4.9 


f 1.9 
-2.4 


2 0+ 1 - 2 


i 7 +3.3 
1 -'-0.8 


1.71};? 



^ This is the light curve shown in Fig. [JJ including both the V 
band data and 5100A continuum AGN Watch data. 



a role in producing the optical variations. To test this hy- 
pothesis, we use the observed X-ray light curve to construct 
a model reprocessed optical light curve, which we can com- 
pare to the observed optical light curve. 

The model, descr ibed in detail by 

iKazanas fc Navakshinl (1200 ll ). assumes that the disc 
temperature is determined by the combined effects of 
viscous dissipation in the disc and variable heating by an 
X-ray source at height h x above the disc and on its axis of 
symmetry. The temperature of the disc is then given by 
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T(R,t)-- 



A)L x {t) 




8naR 2 R 3 



4ira(hi + R 2 ) 



assuming that each annulus of the disc emits as a blackbody. 
Here M is the disc accretion rate in kg s _1 , R the distance of 
the disc surface element from the centre, 7? in the location of 
the inner edge of the accretion disc, L x (t) the variable X-ray 
luminosity, A is the average disc albedo and 9 the angle be- 
tween the disc surface and the direction of the X-ray source. 
R, R ia and h x are measured in units of the gravitational ra- 
dius R g = GM/c 2 , where M is the mass of the black hole in 
kg. The delay between the variations of different wavebands 
is then due to the difference in light travel time to regions 
of the disc of different temperature. 

We show in Figure [9] the lags predicted by this 
model, plotted as a dashed line, assuming M = 1.91 x 
10 6 M Q ((Peterson et al.ll2004l ), M = 0.15M Edd l|Woo fc Urrvl 
120021 ), R„ = 6, h x = 10, A = 0.3 and a constant X-ray lu - 
minosity of L x = 3 x 10 42 erg s" 1 l|Smith fc Vaugharj|2007h . 
The points with error bars are the measured centroid lags 
(of the light curves shown in Figure [2] versus the X-rays) as 
a function of wavelength. With the above fixed parameters, 
the irradiated disc model (Eqn. is not a particularly good 
fit. Although, as we show next, a much better fit can be ob- 
tained if we alter the parameters (e.g. decrease the mass or 
increase the accretion rate). 

Note that both terms in Eqn.[7]scale as R~ 3 ^ for large 
radii. 'Large' here implies R S> -Ri„ and R ^> h x , since we 
can write cos# = h x /(h 2 + R 2 ) 1 ^ 2 . Together with the as- 
sumption that the peak of the emission scales as T oc A -1 , 
and writing the delay between different bands as r = R/c, 
we obtain the well-known r cx A 4 / 3 relationship commonly 
quoted in the context of reprocessing models. (Of course, 
where the large radius assumptions break down the relation- 
ship between the lag a nd the emitte d wavelength is some- 
what modified; see also lGaskertl200sh . We show in Figure 
(dotted line) that a simple A 4 ' 3 model, fitted to the mea- 
sured lags via an arbitrary scaling constant, fits the data 
well. The factor of ~2 difference between the two curves 
can easily be accounted for within the uncertainties on the 
model parameters. 

To construct a model light curve, the emitted flux as a 
function of time may be computed by integrating the emis- 
sion over all light travel time delays r to different parts of 
the disc, and throughout all disc radii, 



hit) 



-Bin 



dR A(R,r)B x (T(R,t-T)) dr 



(8) 



where T(R, t) is given by Eqn. [7] A is an area function ob- 
tained by considering the effect of the disc inclination on 
the distance betwe en the X-ray source and the disc (see 
iBerklev et al"1l2000f ), and B X (T) is the Planck function. R mt 
is a large value, appropriate for the outer edge of the disc. 

We allowed the inner disc radius R ia , the accretion rate 
m = M /Al Edd , height of the X-ray source h x and disc in- 
clination 6 to vary, to find the parameters best describing 
the observed optical light curve. The mass was fixed at the 
rever beration-mapped value, 1.91 x 1O 6 M0 (|Peterson et al.l 
12004 ) . 

Figure [10] shows the observed optical light curve com- 
pared to a typical model reprocessed light curve. Due to the 
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Figure 9. Centroid lags measured from the CCFs between the 
X-rays and optical colours. The dashed line shows the predicted 
lags from an irradiated accretion disc model for parameters ap- 
propriate to NGC 4051 (see text). The dotted line is a simple 
A 4 / 3 model, fitted to the data by an arbitrary scaling constant. 
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Figure 10. Top panel: Galaxy-subtracted V band light curve. 
Middle panel: Best-fit model reprocessed V band light curve, us- 
ing the observed X-ray light curve as input. Bottom panel: Resid- 
ual variations in the optical light curve after the reprocessed com- 
ponent has been removed. 



short light-travel time to the optical emitting region, the 
rapid variability of the X-rays is preserved in the model light 
curves and we find that a range of input parameter values 
produce nearly identical model light curves. As expected, 
the model reprocessed light curve is highly correlated with 
the X-ray light curve, trailing it by 0.4 — 0.8 days, depending 
on the input parameters. 

This is a very simple model, which assumes the X-rays 
to be emitted by a compact source located on the symmetry 
axis of the system, at a height h x above the disc surface. In 
reality, it is more likely that the X-rays come from an ex- 
tended corona. The compact source model used here may 
be viewed as a first order approximation of the centroid 
of such a corona. We believe that this simplified model is 
sufficient in this case, as h x is much smaller than the ra- 
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dius where most of the optical emission is coming from. The 
break in the X-ray power spectrum, which is thought to be 
associated with the inner edge of the disc, corresponds to 
the dynamical time-scale at 35R g , or the viscous time-scale 
of a thick disc (a = 0.1, H/R ~ 0.5) at 3R g , so the bulk 
of the X-rays are probably emitted within a radius smaller 
than, or comparab le to this. Using Comptonization models 
lUttlev et all (2000) estimate the size of the emission region 
to be < 20R g . Integrating the blackbody emission through 
the disc, we find that 95% of the V band emission originates 
from outside a radius of 400-R 9 . From this distance, even an 
extended corona would resemble a point source, making a 
more detailed model of the X-ray source unnecessary. 

Reprocessing by a standard, thin disc appears to de- 
scribe the observed lag and the short time-scale (~days) 
variations of the optical light curve well, but it fails to si- 
multaneously reproduce the larger amplitude flares on time- 
scales of ^months. Closer inspection of the light curves in 
Figure [T] confirms that it is these larger amplitude fluctu- 
ations which lag the X-rays by ~ 40 days (see e.g. MJD 
~ 50750, 52500 and 53700 for obvious examples) . Attempt- 
ing to isolate these variations, we show in the bottom panel 
of Figure 1101 the residual of subtracting the reprocessing 
model light curve from the observed light curve. Cross cor- 
relation of the X-rays with this "residual curve" yields the 
CCF shown in Figure 1111 As the model light curve is very 
well correlated with the X-rays at near days, subtracting 
it removes this correlation peak entirely from the CCF (the 
correlation at this lag drops to zero). Isolating the second 
peak this way, allows us to make a more precise measure- 
ment of the lag it corresponds to. The peak is found to be 
at 38.9I3 3 days, and the 85% centroid at 38.5I2 2 days. 
We calculate the significance of the correlation as before, 
by cross-correlating random red noise light curves with the 
difference light curve. It is found to be significant at 98.4% 
confidence. It is therefore clear that reprocessing by a flat 
disc can only account for the short lags in the light curve and 
that the second peak must be produced by optical emission 
from a structure further out. 



6 DISCUSSION 

The detection of two significant peaks in the X-ray /optical 
correlation function suggest that the optical variability orig- 
inates from more than one location in the system, or through 
more than one process. In the previous section we discussed 
the correlation at short lags in terms of reprocessing by an 
accretion disc. We showed that this model can only repro- 
duce the short time-scale variations in the light curve, and 
that the second peak must be the result of optical emission 
coming from beyond the disc and BLR. 

In this section we investigate the origin of the second 
peak in the correlation function. The lag is found to be con- 
sistent with the light travel time to the inner radius of the 
dust torus in this source, so we investigate the possibility of 
optical emission originating from the torus. Firstly we will 
consider reprocessing of X-rays by the dust, and secondly 
the dust reflecting optical light from the accretion disc. 
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Figure 11. Cross correlation of the reprocessing-subtracted 
residual light curve with the X-rays. The ~ day peak is re- 
moved entirely from the CCF, leaving only the longer time-scale 
correlation peak. The horizontal lines are the 50% (dotted), 95% 
(dot-dashed) and 99% (dashed) significance levels. The peak is 
significant at 98.4%. 

6.1 Second peak in the CCF — optical emission 
from the torus? 

The standard model of AGN holds that the disc and BLR 
is surrounded by a large molecular torus. It is now widely 
accepted that the infrared emission from AGN is the result 
of reproces sing of high energy radiation by hot dust in the 
torus (e.g. I Glass! 120041 ') . We investigate here whether it is 
possible that the second peak detected in the correlation 
function could be the result of optical emission coming from 
the torus. 



6.1.1 Distance to the torus 

The inner radius of the dust distribution is set by the dust 
sublimation radius, inside which the flux from the central 
source will destroy the dust particles. Detailed radiative 
transfer calculations (|Nenkova et al.ll2008h places this radius 
at 

( L \°' 5 / T \~ 26 
r - = °- 4 U 45 ergs-i j (l500Kj PC ' (9) 

where T is the dust temperature and L the UV-optical lu- 
minosity. 

Where adequate data exist, the optical to near-infrared 
delays show reasonable agreement with this model. Us- 
ing both their own monitoring data as well as results 
from l itera ture, Oknvanskii & Hornel l|200lh . lMinezaki et al.l 
( 2004) and lSuganuma et al. (2006) show that the optical-to- 
infrared lag is proportional to L 5 over a range of luminosi- 
ties. 

The exact temperature at which the dust in the torus 
sublimates is not well known, as it depends on the compo- 
sition of the dust. The dust is generally assumed to be a 
mixture of graphite and silicate particles, with sublimation 
temperatures in the range 1500-2000 K. The emission from 
hot dust at the inner edge of the torus is therefore expected 
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to peak in the K band. As graphite can survive higher tem- 
peratures, it is assumed that the dust at the inner edge will 
mos tly cons i st of g raphite grains. 

iGaskelll (|2007l ) pointed out that at such high tempera- 
tures, there must be a small, yet significant amount of optical 
emission coming from the torus as well, as part of the Wien 
tail of the ^1500-2000 K emission. He proposed that the 
increase in time delay measured between the optical bands 
could be the result of contamination by the additional opti- 
cal emission coming from large radii. Optical emission com- 
ing from the torus will of course also be delayed with respect 
to the X-rays, by the same amount as the infrared emission. 
This will cause an asymmetry in the X-ray-optical correla- 
tion function, or, if the signal is strong enough, another peak 
may be seen, corresponding to the light travel time to the 
inner edge of the torus. The dust recombination time-scale is 
expected to be long c ompared to the sublimation time-scale 
l|Koshida et al.l 2009), so that the inner edge of the torus 
is determined by the largest of the X-ray/UV flares. Any 
variations we see in the emitted flux is then probably due 
to temperature variations of the dust (below the sublima- 
tion temperature), rather than a change of the torus inner 
radius. 

The H — K col our measured for NGC 4051 by 
ISuganuma et al.l l|2006l ). corresponds to a blackbody tem- 
perature of approximately 1600 K, so we adopt this value 
as the temperature of the dust in this source. We scale the 
ionizing X-ra y-UV-optical lumi nosity (~13.6 eV-13.6 keV), 
estimated bv lOgle et al.l l|2004) . to a distance of 15.2 Mpc, 
giving L ion = 1.1 x 10 43 erg s _1 . The sublimation radius 
(Eqn. |9j) of NGC 4051 is then r subl = 42.2 light-days. The 
position of the second peak in our correlation function agrees 
well with this value. We note , howe ver, that the V-Xo-K 
lag found by ISuganuma et al.l (2006) is only ~ 20 light- 
days. This may be expected if the disc V band emission 
is contaminated by optical emission from large radii (such 
as the torus), as this will shift the peak of the V vs. K 
CCF to shorter lags. Several studies (|Peterson et al 1 l2000l ; 
iKaspi et al.ll2000l ; iPeterson et al.ll2004l) measure the location 
of the BLR to be at ~ 6 light-days, placing even the smallest 
estimates of the torus well outside the disc and BLR. 

6.1.2 Energetics 

Deta il ed models of AGN dust torii (e.g. iNenkova et alj 
2008a, INenkova et al.l [2008^ 3^ take account of the fact that 
the torus probably consist of clumpy material rather than 
a smooth distribution of dust. One effect the clumps has is 
to cause a steep radial temperature gradient in the torus, as 
optically thick clouds closer to the source can shield those 
further away from direct irradiation. It also takes into ac- 
count the diffuse heating of clumps by each other. Such a 
detailed analysis is beyond the scope of this paper, but as we 
are mainly interested in the hottest dust here, a simplified 
approach may be sufficient to determine whether it is ener- 
getically possible to get the 1-2 mjy V band flux required 
to explain the large amplitude optical fluctuations, from the 
torus. 

A blackbody at 1600 K emits roughly 16% of its total 
emission in the K band, and only 0.04% in the Wien tail 
at V band wavelengths. The fraction of V band emission 
decreases very rapidly as the dust temperature decreases, 



so we do not expect to find a significant contribution to 
the optical from further into the torus. We will assume that 
the optical emission is identified with the inner edge of the 
dust distribution only, where the hottest dust is located. 
For simplicity, we also assume a sharp boundary, although 
the transition is most likely a gradual one, given the differ- 
ence in sublimation temperat ures of diff e rent d ust particles. 
We follow the prescription of iBarvainisI (| 1987f ) to calculate 
the luminosity emitted by the dust. We assume a grain ra- 
dius of a = 0.05/im and infrared (K band) absorption ef- 
ficiency Q a b B = 0.058. The luminosity of a single grain is 
given by Lf, r = A-wa 2 itQ v B„(T) erg s _1 Hz -1 , where B v is 
the Planck spectrum of a grain at temperature T. The K 
band spectral luminosity of a graphite grain at T = 1600 K 
is then = 3.6 x 10~ 17 erg s" 1 Hz " 1 . The average galaxy - 
subtracted K band flux measured bv lSuganuma et al.1 (|2006l ) 
is 31 mjy, which corresponds to a total K band luminosity 
of Lk = 8.6 x 10 27 erg s _1 Hz -1 . Hence the total number 
of hot dust grains is N ~ 2.4 x 10 44 . For UV-optical radia- 
tion, the absorption efficiency is Q aba ~ 1 for grains of this 
size dWickramasinghe et al.l Pl974h . The optical emission of 
a grain at 1600 K is = 1.9 x 10" 19 erg s^ 1 Hz" 1 . If 
we assume the same number of grains to be emitting this 
optical emission (i.e. only the hottest dust) then the to- 
tal V band luminosity we may expect from the torus is 
Lv = 4.4 x 10 25 erg s" 1 Hz -1 . This translates to an ob- 
served flux of only 0.16 mjy, which is an order of magnitude 
smaller than required to explain the additional variability in 
the observed V band light curve. 

One of the main uncertainties in this model is the 
dust composition, as there is no guarantee that dust in 
the vicinity of the active nucleus will be the same as 
dust found in molecular clouds. The temperature of the 
grains is determined by the balance between UV absorp- 
tion and infrared emission, which vary with grain size and 
type. This, in turn, affects the sublimation temperature. 
If we allow the dust temperature to increase to 1800 K 
(which is still within range of the sublimation temperature 
of graphite), the emitted luminosity per grain increases to 
Ly = 1.15 x 10~ 18 erg s _1 Hz -1 , and the additional flux 
observed in the V band to 1.1 mjy. Such an increase in 
the temperature the dust particles can sustain will of course 
move the sublimation radius inwards as well, to a distance 
of r sub i = 31.1 light-days. This is within the la error of the 
peak measured on the X-ray-optical CCF (Figure HU, and 
within 2a of the peak measured from the CCF calculated 
using residual light curve (Figure lll|l . The uncertainty on 
L ion would also allow a factor of ~2 change in the sublima- 
tion radius. A similar analysis for the other optical bands, 
using T = 1800 K and the appropriate absorption efficien- 
cies, yields it = 0.004, B = 0.08, R = 5.8, Rl = 13.2 and 
/ = 16.6 mjy. For wavelengths R and longer, the reprocess- 
ing model therefore predicts a higher flux to come from the 
torus alone than the total galaxy-subtracted flux we mea- 
sure. The measured fluxes in these bands must however in- 
clude a significant contribution from the disc as well, as their 
correlation functions with the X-rays (Fig. [SJ all display a 
significant peak at short lags. 

For increasing wavelength, the torus reprocessing model 
naturally predicts an increase in the amount of reprocessed 
flux originating in the torus, as the redder bands include a 
larger fraction of the total flux emitted by the dust. Our 
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data do not appear to support such an increase, so it is 
unlikely that all the additional optical flux results from dust 
reprocessing. We nevertheless cannot exclude the possibility 
that some of the optical emission in this source is originating 
in the dust torus. 



6.1.3 Optical reflection 

In this section we consider the reflection of optical light by 
the torus. Apart from the illumination by the central source, 
the torus will also receive UV and optical light from the ac- 
cretion disc. The dust absorbs energy at these wavelengths 
very efficiently, so some of this reprocessed optical emission 
will be absorbed by the dust and contribute to the vari- 
able heating of the torus, as discussed in the previous sec- 
tion. Depending on the solid angle subtended by the torus 
at the optically-emitting part of the disc, some of the optical 
emission may also be elastically scattered and reflected by 
the torus. The efficiency of reflection is strongly geometry 
dependent, so we will limit our discussion here to a quali- 
tative description of the main observable results expected. 
Re flection from different geom etries have been considered 
bv lGoosmann fc Gaskelll (2007) in the development of their 
Monte Carlo polarization code STOKES. 

Most of the UV/optical illumination will be relatively 
steady, or will vary on the very long (years) viscous time- 
scales on which the intrinsic emission from the disc will vary. 
However, as we have already shown, a small fraction of the 
optical illumination, resulting from reprocessing of X-rays 
by the accretion disc, will follow the pattern of the X-ray 
variations. As the optically-emitting region is about a light- 
day away from the central X-ray source, the torus will see 
optical emission which broadly follows the X-ray emission, 
but which is smoothed on a time-scale of a day or two rela- 
tive to the X-ray variability pattern. 

The component of the optical emission illuminating the 
torus which has a phase relationship with the X-rays (so 
that it causes a peak in the CCF), must have arisen from 
reprocessing of X-rays in the disc. The amount of scattered 
light depends on the solid angle the torus subtends at the 
optically emitting part of the disc, so for large solid an- 
gles, the elastically scattered optical photons may produce 
a noticeable contribution to the optical variations, and will 
lag th e X-rays by about 40 days. Dust models by iDraind 
(2003) suggest that the scattering albedo (scattering cross 
section as a fraction of the total extinction cross section) of 
Milky Way type dust may be as high as 0.5-0.8, for different 
dust compositions. Dust in AGN probably consists of larger 
grains (smaller grains are evaporated by the intense radia- 
tion field) , so the albedo is uncertain, but it may still be high 
enough for the scattered light to contribute significantly to 
the observed optical radiation. 

Reflection models are supported by spectropolarimet- 
ric studies which are used to detect the broad line region 
in Seyfert 2 galaxies, such as NGC 1068, through polar- 
ized light. In these systems, the continuum source and BLR 
clouds are hidden from our view by the molecular torus, but 
their emission is scattered into our line of sight by warm elec- 
trons and dust lying outside the torus opening, polarizing 
the light. In Seyfert 1 galaxies, however, the optical polariza- 
tion position angle is generally observed to be parallel to the 
radio axis (which is assumed to be the system axis) , suggest- 



ing that the scattering medium has an equatorial geometry, 
perh aps a ring of material co-planar with with accretion 
disc (Smith et al. 2002, 2004), or the inner edge of the torus 
(|Cohen fc Martel2002h . The op tical continuum pola rization 
in NGC 4051, as measured by ISmith et al] l|2002l ). is low 
however. They find 0.5 per cent of the optical continuum at 
7000A to be polarized, amounting to a mere ~ 0.05 mjy. It 
is also not clear to what extent this is intrinsic to the source 
or whether the measurement is contaminated by foreground 
interstellar polarization. If intrinsic, the position angle of the 
polarization is approximately parallel to the radio axis, sug- 
gesting that the scattering occurs in the equatorial region. 
Unfortunately their polarization observation was made dur- 
ing a time for which we do not have optical monitoring data 
available, so we cannot comment on the general optical be- 
haviour of the source at that time. 

There is also evidence for diffuse continuum emission 
and reprocessing by the BLR, as well as reflection by the 
broad line clouds jKorista fc Ferland|[l998l ; lKorista fc Goadl 
l200ll ; lArevalo et al.ll2o'o9f i" in some Seyfert galaxies. Though, 
as mentioned before, any lags associated with the BLR (~6 
days for NGC 4051) are likely to be much shorter than those 
associated with the torus. 

As an aside, we note that the high energy emission 
may also be reflected by the torus. X-ray spectra show ev- 
idence of X-ray reflection by high density gas in the form 
of the "Compton bump" and the Fe Ka fluorescence line 
at 6.4 keV. This line is thought to originate in the accre- 
tion disc, and is therefore expected to be broadened by rel- 
ativistic Doppler motions. Broad Fe lines like this are seen 
in X-ray spectra, but the line is found to have a narrow 
component as well, which is associated with reflection by 
dist ant, optically thick m aterial such as the molecular torus 
(e.g. iNandra et al"1l2007l ). 

Neither the reprocessing nor the scattering model can 
on its own fully account for the optical flux causing the 
second peak in the correlation function, but it is possible 
that either or both may contribute to the optical emission 
observed from the system. Overall we conclude that there 
are many possible sources of optical emission in the inner 
parts of AGN but all those associated with the torus, albeit 
reprocessing or simple reflection, will give more or less the 
same lag. Our current data is not suited to constrain these 
models further, so a more detailed discussion is beyond the 
scope of this paper. 



6.2 The complex optical variability in NGC 4051 

Sh emmer et al.l (|2003l ) find that part of the optical emission 
leads the X-ray variations by approximately 2 days, con- 
trary to what we find here (see Table [3} . As was shown in 
previous sections, there is more than one effect contribut- 
ing to the optical variations in this source. Variations on 
time-scales longer than, or similar to, their 3 month moni- 
toring period may cause asymmetries in the correlation func- 
tion a nd hence a centroid at negative lags. Sh emmer et al.l 
(2003) interpret their result as a combination of reprocess- 
ing and longer time-scale variations, similar to what we find 
over much longer time-scales. We note that the peak of their 
cross-correlation function is consistent with the peak lag we 
find here, ~ 1 day. 

In red noise time series, such as the X-ray and optical 
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light curves used here, peaks and troughs may sometimes 
line up even if the light curves are not really correlated, 
producing spurious peaks in the CCF. We believe that the 
second peak we measure in the CCF here, is not simply due 
to statistical fluctuations like these, as it is present in most 
of the individual segments used to calculate the final CCFs. 
If this was an effect only due to a random correlation in a 
red noise process, we would expect to see it in, e.g. only one 
of the segments' CCFs. We note that the peak does shift be- 
tween ~ 30 — 50 days in the CCFs of the different segments. 
The probability that the peak is due to such statistical fluc- 
tuations was assessed through Monte Carlo simulations in 
Section |4l and was found to be 1.6%. 

So far we have interpreted the second peak in the CCF 
only in terms of a light crossing time-scale. Other time-scales 
of interest are the dynamical (Kepler) time-scale, and the 
viscous time-scale, on which we expect accretion rate fluc- 
tuations to propagate through the disc. A lag of the optical 
variations behind the X-rays, as seen here, would require 
such fluctuations to be propagating outwards through the 
disc for this time-scale to be relevant. A viscous time-scale of 
40 days corresponds to R = 50R g in a moderately thin, op- 
tically thick disc (H/R ~ 0.1, a di3C = 0.1), or to R = 4407? 9 , 
in a geometrically thick disc (H/R ~ 0.5, Q d i S c = 0.1). 
R = 50R g lies well within the Far- to Extreme UV emit- 
ting part of the disc; less than 0.5% of the optical emission 
come from within this radius. R = 440i? g contains approxi- 
mately 6% of the optical emission. Enhanced emission due to 
viscous instability in this region could therefore only cause 
a ~ 6% variation in the light curve. This translates to a 
flux variation of 0.2 mjy, which is not enough to explain the 
variations seen in the light curve. As an orbital time-scale, 
40 days corresponds to the radius which contains 69% of 
the V band emission, but is is difficult to associate an or- 
bital time-scale with a lag between the emission in different 
wavebands, as we see here. 

If this is indeed optical emission from the torus, why 
has this not been seen in other sources? NGC 4051 is a low 
mass, low luminosity system, which means that it is physi- 
cally smaller than higher mass AGN. Lags between the op- 
tical and near-infrared suggest that dust sublimation radii 
in most sources are several tens to hundreds of days, but 
here it is only ~20-40 days. The optical emission from the 
torus will be a very weak signal, so long enough light curves 
are necessary to include several cycles of the correlated vari- 
ations to produce a significant peak in the CCF. We note 
that calculating the CCF using only half of the V band data, 
decreases the significance of the second peak in Figure [6] to 
96.3%, or to below 95% if only a quarter of the light curve 
is used. 

More importantly, in low mass systems the V band 
emission originating in the torus is probably a larger frac- 
tion of the total V band emission emitted by the system, 
than it is in more massive systems. This is due to the varia- 
tion in disc temperature with black hole mass. For example, 
for a fixed accreti on rate of rh = 0.10, a standard disc (e.g. 
iFrank et all 120021 ') around a 10 8 M Q black hole emits 1.9% 
of its total luminosity at V band wavelengths, but for a 
1O 6 M0 black hole the disc is hotter, the bulk of the energy 
is emitted at bluer wavelengths and the fraction of the total 
emission emitted at V band wavelengths decreases to 0.4%. 
On the contrary, the torus in all sources is expected to be at 



approximately the same temperature (as its inner radius and 
temperature is determined by dust sublimation) and hence 
it will contribute a constant fraction of the V band emission. 
Comparing this constant fraction with the V band emission 
originating in the disc, the torus contributes a larger frac- 
tion to the total V band emission in the 10 6 Mq source than 
it does in the 1O 8 M0 source. In more massive systems the 
V band emission is therefore dominated by the disc, so the 
weak optical emission from the torus will be very hard to 
isolate. 

Long-term multicolour optical and infrared monitoring 
of NGC 4051 is needed to reliably measure the inner radius 
of the dust torus and t he possible contribu tion of the disc to 
the near-infrared (e.g. iTomita et al1l2006h . Detailed model- 
ing of the torus emission, taking into account both the ge- 
ometry of the system as well as the dust albedo and emission 
properties, is also necessary to properly assess the contribu- 
tion this could have to the optical emission observed. 



7 SUMMARY 

We have presented a 12.6 year 3-10 keV X-ray light curve 
from our ongoing long term X-ray monitoring program 
on RXTE, along with concurrent optical monitoring light 
curves. We combined the V band data from four different 
telescopes with archival 5100A continuum data, to construct 
an optical light curve spanning the same length in time. We 
also include light curves in the u, B, R, Rl and / optical 
bands, each of length between 3.6 and 6.6 years in time. 

We find both the X-ray and optical emission to be 
highly variable, but find the amplitude of the optical emis- 
sion to be smaller than that of the X-rays, even after sub- 
tracting the contaminating host galaxy flux falling inside 
the aperture. The power spectra show that the optical vari- 
ability power is lower than the X-rays on all time-scales 
probed. With the current data, the optical PSD is best de- 
scribed by an unbroken power law of slope a — 1.4, which 
is much steeper than the low frequency slope of the X-ray 
PSD (oil = 1.1). This suggests that the variations in the 
different bands do not have a common origin or that they 
are not related in a simple way. We will address this issue 
in more detail in a future paper. 

We measured the galaxy flux inside the 15" aperture by 
decomposition of the optical images, and find the measured 
fluxes to agree well with a template spectrum of a similar 
non- active spiral galaxy. 

The cross-correlation function between the optical light 
curve and the X-rays displays two statistically significant 
peaks. The strongest peak corresponds to the optical varia- 
tions lagging the X-rays by 1.2lJ'g days. This time-scale is 
consistent with the light travel time to the optical emitting 
region of the accretion disc, and points to X-ray reprocess- 
ing as the source of the fast optical variations. The delays 
measured between other optical colours and the X-rays are 
also consistent with this interpretation. However, using a 
reprocessing model, we find that all the optical variability 
cannot be accounted for this way. Reprocessing appears to 
describe the small amplitude, fast variations well, but it fails 
to reproduce the larger amplitude variations on time-scales 
of months-years that are also seen in the light curve. 

The second significant peak in the CCF is at an opti- 
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cal lag of 38.9tg 4 days behind the X-rays. By subtracting 
a model reprocessed light curve from the observed optical 
light curve and recalculating the cross-correlation, we can 
constrain this peak more precisely, viz. 38.9I33 days. Inter- 
preting this as a light travel time, the corresponding radius 
is consistent with the dust sublimation radius in this AGN, 
which is a theoretical estimate of the inner radius of the dust 
torus. 

As such, we have investigated the origin of the addi- 
tional optical emission, both in terms of a weak optical con- 
tribution from X-ray reprocessing by the dust, and in terms 
of reflection of optical light from the disc by the dust. We 
find that neither process can on its own explain the origin 
of the second peak in the CCF of NGC 4051. It remains 
possible however, that both these processes may contribute 
to the optical emission observed. Detailed modeling of the 
dust emission properties at optical wavelengths, as well as 
further optical polarimetry will be required to thoroughly 
test this hypothesis. 

In general, the weak (though significant) short time- 
scale correlation reported here is consistent with a model 
where the strength of the correlation is determined by 
the mass of the ce ntral black hole and disc accretion rate 
jUttlev et al.ll2003l ). When scaled in te rms of gravitational 
radiu s, a standard accretion disc (e.g. IShakura fc Sunvaevl 
Il973h around a smaller mass black hole, with a given accre- 
tion rate, is hotter. The optical emitting region will therefore 
be further away from the centre and the viscous time scale 
here will be much longer than near the centre, where the 
X-rays are emitted. The variations in the different bands 
may therefore be expected to vary incoherently, weaken- 
ing the strength of the correlation. Furthermore, the mass- 
dependence of the disc temperature may also explain why an 
optical contribution from the torus is detected in this source 
and not in more massive systems: the cooler discs around 
more massive systems will emit most of their luminosity at 
optical wavelengths, dominating over the weak optical sig- 
nal from the torus. In low mass systems with hotter discs, 
such as NGC 4051, the peak of the disc emission is shifted 
to shorter wavelengths, decreasing the fraction of the emis- 
sion emitted at V band wavelengths and hence increasing 
the possibility of detecting optical emission from the torus. 



ACKNOWLEDGEMENTS 

We thank the referee for constructive suggestions which im- 
proved this paper. 

We are grateful to the staff of the Rossi X-ray Timing 
Explorer, the Liverpool Telescope and the Faulkes Telescope 
for their support of our long term program and also thank 
Keith Home and Ed Cackett for their contribution in ob- 
taining this data. 

EB gratefully acknowledges financial support in the 
form of a Stobie-SALT Scholarship from South African 
National Research Foundation and the University of 
Southampton. IM C H acknowledge support from the STFC 
under rolling grant PP/D001013/1 and PA acknowledges 
support from the Natural Science Foundation of China 
(grants 10773024, 10833002, 10821302, 10825314, 10873030 
and 40636031), and National Basic Research Program of 
China (grant 2009CB824800). PU acknowledges support 



from an STFC Advanced Fellowship. The work of SGS and 
NGC is partially supported by the CosmoMicroPhysics Tar- 
get Complex Program of the Ukrainian NAS. 



REFERENCES 

Arevalo P., Uttley P., 2006, MNRAS, 367, 801 

Arevalo P., Uttley P., Kaspi S., Breedt E., Lira P., 

McHardy I. M., 2008, MNRAS, 389, 1479 
Arevalo P., Uttley P., Lira P., Breedt E., JVTHardy I. M., 

Churazov E., 2009, MNRAS, 397, 2004 
Barvainis R., 1987, ApJ, 320, 537 

Bentz M. C, Peterson B. M., Netzer H., Pogge R. W., 

Vestergaard M., 2009, ApJ, 697, 160 
Berkley A. J., Kazanas D., Ozik J., 2000, ApJ, 535, 712 
Bica E., 1988, A&A, 195, 76 

Breedt E., Arevalo P., McHardy I. M., Uttley P., Sergeev 
S. G., Minezaki T., Yoshii Y., Gaskell C. M., Cackett 
E. M., Home K., Koshida S., 2009, MNRAS, 394, 427 

Cackett E. M., Home K., Winkler H., 2007, MNRAS, 380, 
669 

Chesnok N. G., Sergeev S. G., Vavilova I. B., 2009, Kine- 
matics and Physics of Celestial Bodies, 25, 107 

Cohen M. H., Martel A. R., 2002, in Crenshaw D. M., 
Kraemer S. B., George I. M., eds, Mass Outflow in Active 
Galactic Nuclei: New Perspectives Vol. 255 of Astronom- 
ical Society of the Pacific Conference Series, Spectropo- 
larimetry of Broad Ha Lines and Geometry of the BLR. 
p. 255 

Collier S. J., Home K., Kaspi S., Netzer H., Peterson 
B. M., Wanders I., Alexander T., Bertram R., Comas- 
tri A., Gaskell C. M., Malkov Y. F., Maoz D., Mignoli 
M., Pogge R. W., Pronik V. I., et al 1998, ApJ, 500, 162 

Done C, Ward M. J., Fabian A. C, Kunieda H., Tsuruta 
S., Lawrence A., Smith M. G., Wamsteker W., 1990, MN- 
RAS, 243, 713 

Doroshenko V. T., Sergeev S. G., Merkulova N. I., Sergeeva 
E. A., Golubinsky Y. V., Pronik V. I., Okhmat N. N., 
2005, Astrophysics, 48, 156 

Draine B. T., 2003, ApJ, 598, 1017 

Frank J., King A., Raine D. J., 2002, Accretion Power in 
Astrophysics: Third Edition. Accretion Power in Astro- 
physics, by Juhan Frank and Andrew King and Derek 
Raine, pp. 398. ISBN 0521620538. Cambridge, UK: Cam- 
bridge University Press, February 2002. 

Gaskell C. M., 2007, in Ho L. C, Wang J.-W., eds, The 
Central Engine of Active Galactic Nuclei Vol. 373 of As- 
tronomical Society of the Pacific Conference Series, The 
Origin of Wavelength-dependent Continuum Delays in 
AGNs - A New Model, p. 596 

Gaskell C. M., 2008, in Revista Mexicana de Astronomia y 
Astrofisica Conference Series Vol. 32 of Revista Mexicana 
de Astronomia y Astrofisica Conference Series, Accretion 
Disks and the Nature and Origin of AGN Continuum Vari- 
ability, pp 1-11 

Gaskell C. M., Sparke L. S., 1986, ApJ, 305, 175 

Glass I. S., 2004, MNRAS, 350, 1049 

Goosmann R. W., Gaskell C. M., 2007, A&A, 465, 129 

Haardt F., Maraschi L., 1991, ApJL, 380, L51 

Kaspi S., Smith P. S., Netzer H., Maoz D., Jannuzi B. T., 
Giveon U., 2000, ApJ, 533, 631 



X-ray/optical variability in NGC 4051 17 



Kazanas D., Nayakshin S., 2001, ApJ, 550, 655 
Klimek E. S., Gaskell C. M., Hedrick C. H., 2004, ApJ, 
609, 69 

Koratkar A., Blaes O., 1999, PASP, 111, 1 

Korista K., Ferland G., 1998, ApJ, 495, 672 

Korista K. T., Goad M. R., 2001, ApJ, 553, 695 

Koshida S., Yoshii Y., Kobayashi Y., Minezaki T., Sakata 

Y., Sugawara S., Enya K., Suganuma M., Tomita H., Aoki 

T., Peterson B. A., 2009, ApJL, 700, L109 
Krolik J. H., Home K., Kallman T. R., Malkan M. A., 

Edelson R. A., Kriss G. A., 1991, ApJ, 371, 541 
Lyubarskii Y. E., 1997, MNRAS, 292, 679 
Maoz D., Edelson R., Nandra K., 2000, AJ, 119, 119 
Maoz D., Markowitz A., Edelson R., Nandra K., 2002, AJ, 

124, 1988 

Mason K. O., McHardy I. M., Page M. J., Uttley 
P., Cordova F. A., Maraschi L., Priedhorsky W. C, 
Puchnarewicz E. M., Sasseen T., 2002, ApJL, 580, L117 

McHardy I. M., Koerding E., Knigge C, Uttley P., Fender 
R. P., 2006, Nature, 444, 730 

McHardy I. M., Papadakis I. E., Uttley P., Page M. J., 
Mason K. O., 2004, MNRAS, 348, 783 

Minezaki T., Yoshii Y., Kobayashi Y., Enya K., Suganuma 
M., Tomita H., Aoki T., Peterson B. A., 2004, ApJL, 600, 
L35 

Nandra K., Le T., George I. M., Edelson R. A., Mushotzky 
R. F., Peterson B. M., Turner T. J., 2000, ApJ, 544, 734 

Nandra K., O'Neill P. M., George I. M., Reeves J. N., 2007, 
MNRAS, 382, 194 

Nenkova M., Sirocky M. M., Ivezic Z., Elitzur M., 2008, 
ApJ, 685, 147 

Nenkova M., Sirocky M. M., Nikutta R., Ivezic Z., Elitzur 

M., 2008, ApJ, 685, 160 
Ogle P. M., Mason K. O., Page M. J., Salvi N. J., Cordova 

F. A., McHardy I. M., Priedhorsky W. C, 2004, ApJ, 606, 

151 

Oknyanskij V. L., Home K., 2001, in Peterson B. M., Pogge 
R. W., Polidan R. S., eds, Probing the Physics of Ac- 
tive Galactic Nuclei Vol. 224 of Astronomical Society of 
the Pacific Conference Series, Reverberation Radii of Dust 
Holes in Active Galactic Nuclei, p. 149 

Oknyanskij V. L., Home K., Lyuty V. M., Sadakane K., 
Honda S., Tanabe S., 2003, in Collin S., Combes F., Shlos- 
man I., eds, Active Galactic Nuclei: From Central Engine 
to Host Galaxy Vol. 290 of Astronomical Society of the Pa- 
cific Conference Series, Time Delays in NGC 4151. p. 119 

Peng C. Y., Ho L. C, Impey C. D, Rix H.-W., 2002, AJ, 
124, 266 

Peterson B. M., 1999, in Gaskell C. M., Brandt W. N., Diet- 
rich M., Dultzin-Hacyan D., Eracleous M., eds, Structure 
and Kinematics of Quasar Broad Line Regions Vol. 175 
of Astronomical Society of the Pacific Conference Series, 
The International AGN Watch: Past, Present, and Future, 
p. 49 

Peterson B. M., Ferrarese L., Gilbert K. M., Kaspi S., 
Malkan M. A., Maoz D., Merritt D., Netzer H., Onken 
C. A., Pogge R. W., Vestergaard M., Wandel A., 2004, 
ApJ, 613, 682 

Peterson B. M., McHardy I. M., Wilkes B. J., Berlind 
P., Bertram R., Calkins M., Collier S. J., Huchra J. P., 
Mathur S., Papadakis I., Peters J., Pogge R. W., Romano 
P., Tokarz S., Uttley P., Vestergaard M., Wagner R. M., 



2000, ApJ, 542, 161 
Russell D. G., 2003, ArXiv Astrophysics e-prints 
Salvati M., Hunt L. K., Calamai G., Del Zanna G., Gi- 

annuzzo E., Kidger M., Mannucci F., Stanga R. M., 

Wamsteker W., 1993, A&A, 274, 174 
Santos J. F. C. J., Alloin D., Bica E., Bonatto C, 2002, in 

Geisler D. P., Grebel E. K., Minniti D., eds, Extragalac- 

tic Star Clusters Vol. 207 of IAU Symposium, A spectral 

library of galaxies, star clusters and stars., p. 1 
Sergeev S. G., Doroshenko V. T., Golubinskiy Y. V., 

Merkulova N. I., Sergeeva E. A., 2005, ApJ, 622, 129 
Shakura N. I., Sunyaev R. A., 1973, A&A, 24, 337 
Shemmer O., Uttley P., Netzer H., McHardy I. M., 2003, 

MNRAS, 343, 1341 
Smith J. E., Robinson A., Alexander D. M., Young S., Axon 

D. J., Corbett E. A., 2004, MNRAS, 350, 140 
Smith J. E., Young S., Robinson A., Corbett E. A., Gian- 

nuzzo M. E., Axon D. J., Hough J. H., 2002, MNRAS, 

335, 773 

Smith R., Vaughan S., 2007, MNRAS, 375, 1479 
Steele I. A., Smith R. J., Rees P. C, Baker I. P., Bates S., 
Bode M., et al. 2004, Ground-based Telescopes. Edited 
by Oschmann, Jacobus M., Jr. Proceedings of the SPIE, 
5489, 679 

Suganuma M., Yoshii Y., Kobayashi Y., Minezaki T., Enya 
K., Tomita H., Aoki T., Koshida S., Peterson B. A., 2006, 
ApJ, 639, 46 

Summons D., 2008, Ph.D. Thesis, University of Southamp- 
ton 

Timmer J., Koenig M., 1995, A&A, 300, 707 

Tomita H., Yoshii Y., Kobayashi Y., Minezaki T., Enya K., 

Suganuma M., Aoki T., Koshida S., Yamauchi M., 2006, 

ApJL, 652, L13 
Uttley P., Edelson R., McHardy I. M., Peterson B. M., 

Markowitz A., 2003, ApJL, 584, L53 
Uttley P., McHardy I. M., Papadakis I. E., 2002, MNRAS, 

332, 231 

Uttley P., McHardy I. M., Papadakis I. E., Cagnoni I., 

Fruscione A., 2000, MNRAS, 312, 880 
Vaughan S., Edelson R., Warwick R. S., Uttley P., 2003, 

MNRAS, 345, 1271 
Wanders I., Peterson B. M., Alloin D., Ayres T. R., Clavel 

J., Crenshaw D. M., Home K., Kriss G. A., Krolik J. H., 

1997, ApJS, 113, 69 
White R. J., Peterson B. M., 1994, PASP, 106, 879 
Wickramasinghe N. C, Lukes T., Dempsey M. J., 1974, 

ApSS, 30, 315 
Woo J.-H., Urry C. M., 2002, ApJ, 579, 530 
Yoshii Y., Kobayashi Y., Minezaki T., 2003, in Bulletin 

of the American Astronomical Society Vol. 35 of Bulletin 

of the American Astronomical Society, The MAGNUM 

(Multicolor Active Galactic NUclei Monitoring) Project. 

p. 752 

This paper has been typeset from a Tp^X/ DTffjX file prepared 
by the author. 



